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ABSTRACT 

We present a detailed study of the physical properties of the molecular gas in a sample 
of 18 molecular gas-rich early-type galaxies (ETGs) from the ATLAS'^^ sample. Our 
goal is to better understand the star formation processes occurring in those galaxies, 
starting here with the dense star-forming gas. We use existing integrated ^^C0(1- 
0, 2-1), 13CO(1-0, 2-1), HCN(l-O) and HCO+(1-0) observations and new ^^CO{3-2) 
single-dish data. From these, we derive for the first time the average kinetic tempera- 
ture, H2 volume density and column density of the emitting gas in a significant sample 
of ETGs, this using a non-local thermodynamical equilibrium (non-LTE) theoretical 
model. Since the CO lines trace different physical conditions than of those the HCN 
and HCO"*" lines, the two sets of lines are treated separately. For most of the molec- 
ular gas-rich ETGs studied here, the CO transitions can be reproduced with kinetic 
temperatures of 10-20 K, H2 volume densities of 10'^"'* cm"'^ and CO column densities 
of 10^^"^° cm~^. The physical conditions corresponding to the HCN and HCO"'" gas 
component have large uncertainties and must be considered as indicative only. We 
also compare for the first time the predicted CO spectral line energy distributions 
(SLEDs) and gas properties of our molecular gas-rich ETGs with those of a sample of 
nearby well-studied disc galaxies. The gas excitation conditions in 13 of our 18 ETGs 
appear analogous to those in the centre of the Milky Way, hence the star formation 
activity driving these conditions is likely of a similar strength and nature. Such results 
have never been obtained before for ETGs and open a new window to explore further 
star-formation processes in the Universe. The conclusions drawn should nevertheless 
be considered carefully, as they are based on a limited number of observations and on 
a simple model. In the near future, with higher CO transition observations, it should 
be possible to better identify the various gas components present in ETGs, as well as 
more precisely determine their associated physical conditions. To achieve these goals, 
we show here from our theoretical study, that mid- J CO lines (such as the CO (6-5) 
line) are particularly useful. 

Key words: Submillimetre, ISM: molecules, galaxies: elliptical and lenticular-cD, 
stars: formation, methods: data analysis 



1 INTRODUCTION 

Early-type galaxies (ETGs) are central to our understand- 
ing of the mass assem bly history of galaxie s. Historically, 
the SAURON projeclQ l|de Zeeuw et al.ll2002l ') aimed to un- 
derstand the formation and evolution of elliptical galaxies, 
lenticular galaxies and spiral bulges from three-dimensional 
(3D) observations, this usin g the panoramic op tical integral- 
field spectrograph SAURON (jBacon et al.l 120011 ). However, a 
more statistically robust sample was required to test its con- 
clusions, and more recently the ATLAS^'^ projecfl has tar- 
geted a complete sample of 260 early-type galaxies within 
the local (42 Mpc) volume, also combining multi-wavelength 
surveys, numerical si mulations and semi-anal ytic modelling 
of galaxy formation ijCappellari et al1l2011al . hereafter Pa- 
per I). ATLAS=^° thus aims to constrain the dynamics of a 
complete volume-limited sample of ETGs, to fully charac- 
terise this class of galaxies and explore their formation and 
evolution. 

The ATLAS™ dataset now provides an accurate in- 
ventory of the baryon budget of local ETGs, with a de- 
tailed two-dimensional (2D) description of the stellar and 
gaseous (both cold and warm) kinematics, and resolved 
stellar population information for all galaxies. It allows to 
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probe the fossil record and ask a number of specific ques- 
tions. How do slow rotators form? What are the physi- 
cal processes determining their kinematic and photomet- 
ric properties? What is the role of major and minor merg- 
ers in their formation history? How is star formation in 
ETGs quenched? Is the mechanism different for fast- and 
slow-rotating ETGs? How does it depend on environment? 
Etc. Earlier ATLAS^^ papers have already partiall y ad- 
dress ed many of these qu estions (e.g. IBoIs et al.i 2010l . 
I2OIII , hereafter P aper VI; ICappellari et all boilU here- 
after Paper VII; Davis et al.l 2011a , here after Paper X; 



, , __. .Davis et al.l . 

Emsellem et al.ll201ll , here after Paper HI; iKhochfar et al.l 
201ll . hereafter Paper VIII; iKrainovic et al.l 20 111 , hereafter 



Paper II; You ng et al.ll201ll . hereafter Paper IV). 

effort to study the cold 



ATLAS=*° has a 
gas in ETGs, both 



built-in 

atomic 

I2OO6I ; lOosterloo et"al] I2OI0I ; 
after Paper XIII) and 



(e.g. HI; 


Morganti et al.l 


ISerra et al. 


2OI2I. here- 



molecular (e.g. 



'CO; 



ICombes. Young fc Bureaul2007l ; Paper IV). These and other 



efforts in the literature (e.g. 


Knapp. Turner & Cunniffe 


1985 


; Kennev & Yound 


19861; Knapp & Rupen 


1996 


; Welch & Sage' '20031; ISage & WelchI |2006|; 


Sage. Welch & Young 2007) have been fuelled bv the 



realisation over the past few decades that many ETGs, 
previously thought to be 'red and dead', actually harbour 
a substantial amount of cold (hence molecular) gas, and 
that many have a non-negligible level of star formation 
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(e.g. I Welch fc Sagd l2003l : ICombes. Young fc Bureaul l2007l : 
ICrocker et al.ll201ll ). 

More chemically complex species such as HCO"*" and 
HCN as well as the less abundant isotopomer ^'^CO have 
also recently been observe d in ATLAS^° ETGs (|Krips et all 
I2OIOI : ICrocker et al.|[20ll . hereafter Paper XI), mostly us- 
ing single-dish telescopes. Paper XI presents IRAM 30m 

telescope observations of the 18 ^^CO-brightest ATLAS^'^ 2.1 Literature data 
ETGs in the ^=^00(1-0), "C0(2-l), HCN(l-O) and HCO+(l- 
0) lines. Several key molecular line ratios are investigated as 
empirical indicators of the molecular gas physical conditions. 
Comparison with spiral galaxies shows that the line ratios of 
ETGs and spirals generally overlap, but many ETG outliers 
also exist. A few ETGs have more starburst-like line ratios, 
while a few others have very low i^CO/^CO ratios. Cor- 
relations of line ratios involving ^^CO(l-O) are found with 
many galaxy properties: molecular-to-atomic gas mass frac- 
tion, dust temperature, dust morphology and stellar popu- 
lation age. Overall, it seems as if the driver of these correla- 
tions is a reduced optical depth in ETGs that have recently 
acquired or are in the process of acquiring their cold gas. 

While these results are suggestive. Paper XI does not 
include modelling of the line ratios to derive the physi- 
cal conditions of the molecular gas. Detailed modelling has 
been performed for only one ETG to date, Centaurus A 



briefly. We remind the reader that all the basic properties 
of our sample galaxies can be found in previous ATLAS^'' 
papers (e.g. Paper I). 



2 OBSERVATIONS 



Most of the i^CO, ^^CO, HCN and HCO+ molecular gas 
detections in ETGs that we are using here for modelling are 
taken from previously published SAURON and ATLAS^° 
follow-ups. Observational parameters b eyond those detailed 
below can be found in these studies (jWelch fc Sagell2003l : 
ICombes. Young fc Bureaul l2007l : iKrips et al.l l2010l : ; Pa- 
per IV; Paper XI). To be more precise, observations of the 
HCN(l-O) {u = 89.087 GHz), HCO+(1-0) {u = 89.188 GHz), 
^^CO(l-O) {u = 110.201 GHz) and "CO(2-l) (u = 
220.398 GHz) lines were presented in iKrips et al.l (|2010l ) 
and Paper XI, following-up on the 18 strongest ^■^CO(l-O) 
{u = 115.2 71 GHz) and ^^CO(2-l) (v ~ 23 0.538 GHz) de- 
tections of lCombes. Young fc Bureaul |20o3) and Paper IV. 



Deen periormea tor only one ij^iu to gate, uentaurus A 
l|Ebneter fc Balick! 11983 : lEckart et al] ll99oi : iLiszj l200j ). 
with several m olecular absorption lines observed towards 
its nucleus fe.g. [Eckart et al.lll990l: jWild. Eckart fc Wikhndl 



ll997l : lMuller fc Dinh-V-Trundl2009l '). No svstematic and ho- 
mogenous study, based on a comparison of observations of 
a well-defined sample of ETGs with theoretical models, has 
ever been performed. To better understand the evolution of 
the cold gas and associated star formation in ETGs, it is 
however essential to constrain gas physical parameters such 
as the kinetic temperature, H2 volume density and the col- 
umn density of species such as CO, HCN and HCO"'', even 
if they are averaged over entire galaxies. 

Here, we therefore aim to estimate the globally-averaged 
physical properties of the molecular gas in a significant sam- 
ple of ETGs (18), using a single-dish dataset (largely the 
observations presented in Paper XI) and focussing on the 
two gas components traced by the CO and HCN/HCO''' 
lines, respectively. As this first study is by necessity ex- 
ploratory, we use a simple theoretical framework, non-local 
thermodynamical equilibrium (non-LTE) model, specifically 
the Large Velocit y Gra dient (LVG) model developed by 
Ivan der Tak et al] |20o3). We have also recently obtained 
^■^ CO (3-2) observations with the HARP instrument on the 
James Clerk Maxwell Telescope (JCMT) for a few objects. 
We report both on these new observations and the theoret- 
ical model in this paper. 

The structure of the paper is as follows. We first present 
new CO (3-2) observations of 8 sample galaxies with existing 
multi-line observations in Section [S] and briefly review the 
main characteristics of the other observations on which our 
theoretical analysis is based. The main assumptions and pa- 
rameters of the LVG model are described in Section O The 
comparison of the model predictions and observations, and 
the identification of the models reproducing best the data 
are detailed in Section U A discussion of our results and 
literature data are then presented in Section [5] The results 
are discussed further in Section [6] where we also conclude 



2.2 New ^^CO(3-2) observations 

New ^^CO(3-2) {u = 345.795 GHz) observations of 8 of 
those 18 galaxies were acquired at the JCMT (project 
mllaull) between 31 March and 4 September 2011. These 
galaxies are pathfinders for a larger ^^CO(3-2) survey, and 
were selected because of their bright ^■^CO(l-O) fluxes. For 
5 galaxies (IC0676, IC1024, NGC6014, PGC058114 and 
UGC09519), the HARP receiver and ACSIS backend were 
used in a conflguration yielding a 1 GHz total bandwidth 
and 2048 channels each 488 kHz-wide. Three galaxies 
(NGC3665, NGC4526 and NGC5866) have broader lines 
(see Fig. [TJ and for those the backend conflguration used 
two overlapping sub-bands, each 1 GHz wide and covering 
a total bandwidth of 1.9 GHz with 977 kHz-wide channels. 
Observations were carried out in the grid-chop or staring 
mode, with beam switching and a 90" throw. Opacities 
at 225 GHz {T225) ranged from 0.05 to 0.1 and system 
temperatures (at 345 GHz) from 230 to 370 K. Pointing 
tests before each galaxy observation and less frequent focus 
calibrations were made on the standard sources RXBoo, 
IRC+10216, IRC+20326, CIT6, RTVir, and XHer, and the 
average pointing offset was 3"-4". Repeated observations 
of the standard sources G34.3 and 16293-2422 every night 
suggest that the absolute calibration is constant at the 
2 to 5% level. Stable baselines required only a O"' order 
baseline subtraction. Conversion from antenna tempera- 
ture to main beam brightness temperat ure was done via 
3mb = TX/0.6 (see 1 Warren e t al. 2 010l or the references 
http://www.jach.hawaii.edu/ JCMT/spe ctral_line/| and 



ihttp: / /www .jach.hawaii.edu/JCMT/spectraLline/Standards' 
/current_cals.html). The half-power beam width (HPBW) 
of the JCMT at uRi 345 GHz is 14". 

The data pre- reduction for the new CO (3-2) observa- 
tions was done using Starlink software (KAPPA, SMURF 
and STLCONVERT packagesEI), and the data were subse- 



^ See |http: / /www.jach.ha waii.edu/ JCMT/spectral.line/l 

dat a_reduction / acsisdr/basics .html 
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Table 1. Velocity-integrated line intensities of the new CO (3-2) 
JCMT data (see Section O. 



Galaxy 






(K km s-i) 


IC0676 


10.14±0.34 


IC1024 


9.57±0.74 


NGC3665 


12.32±0.41 


NGC4526 


18.58±0.58 


NGC5866 


23.84±0.59 


NGC6014 


8.60±0.46 


PGC058114 


9.01±0.47 


UGC09519 


11.93±0.52 



quently translated to CLAS^ format for final reduction. 
The reduced spectra for all 8 galaxies are presented in Fig- 
ure [1] and the integrated intensities are listed in Table [T] No 
other CO(3-2) observation is available so far. The ^^CO(3- 
2) line profile in NGC4526 is asymmetric, probably caused 
by a combination of a small pointing error and the emission 
region being slightly larger than the beam at this frequency 
(causing some "missing" flux at low velocities). With only 
single-dish information, it is difficult to precisely estimate 
how much fiux is missed. However, we do have recently ac- 
quired interferometric ^^00(2-1) CARMA data of NGC4526 
(Davis et al. in prep.), and estimate ~ 8% of missing fiux. 
We have checked that our conclusions remain unchanged 
when using the corrected flux (see Section|4|. The CO (3-2) 
line emission of IC1024 suffers from a greater pointing un- 
certainty, that does not allow us to use this detection in our 
modelling work (Section [3]). Observations of IC1024 alone 
used the jiggle-chop mode with a 4 x 4 map-centered grid, 
yielding a square grid of spectra spaced by 7'.' 5. The optical 
center of the galaxy therefore does not fall on a grid posi- 
tion, and we averaged the four central spectra to produce 
the spectrum shown in Figure [T] The resulting averaged line 
proflle is however signiflcantly different from those of the 
CO (1-0) and CO (2-1) lines (see Paper XI), with a peak 
flux velocity more than 100 km s~^ away from those of the 
lower-J CO fines and a ^^CO{3-2) fine width 20% smaller. 
Any CO (3-2) flux correction would thus be highly uncer- 
tain, and we do not include this detection in our modelling 
work. Overall, we therefore have reliable ^^CO(3-2) obser- 
vations for 7 of our 18 sample galaxies. 



2.3 Data reduction and analysis 

The data reduction for all galaxies, including the ^^CO(3- 
2) line, has been performed (or redone when already pub- 
lished) in a systematic and homogeneous way using the 
GILDAS package, as described in Paper XI. As the line 
profiles are not always Gaussian, velocity-integrated line in- 
tensities were extracted by simply summing the intensities 
over the range of velocities with detected ^'^ CO (1-0) emis- 
sion (see Paper IV). Although the ^^00(3-2) emission of 
NGC6014 and PGC058114 is slightly off-centred within that 



range, we have kept it for consistency with previously pub- 
lished results. The statistical uncertainties resulting from 
this method are discussed in great detail in Paper XI and 
lead to errors on the velocity-integrated line intensities of 
^5 10%. If the peak fiux is lower than three times the rms 
noise in line-free channels (integrated over the expected ve- 
locity range), as is the case for several galaxies in the HCN(1- 
0) and HCO^(l-O) lines, the line is considered undetected 
and an upper limit of three times the rms noise is adopted. 

The LVG models (see Section [3.2|) also require the ob- 
served line widths of each transition and molecule. Depend- 
ing on the source and spectrum, we use either a Gaussian 
or a double-peaked line profile fit to extract the line width, 
taking the line width as the full-width at half maximum 
(FWHM) for Gaussian profiles and the velocity width given 
by the shell method for double-peaked profiles Q These pro- 
vide satisfactory estimates given the crude model assump- 
tions. The line widths are tabulated in Table [2] where the 
errors reported are formal fitting uncertainties, with an av- 
erage of 20%. For undetected lines (flux upper limits), we 
take a line width equal to that of the 'closest' molecular 
species. For example, when undetected, we have assumed 
that HCO+(1-0) has the same line width as the HCN(l-O), 
since we consider that both species trace the same gas com- 
ponent. Similarly, when the ^^CO(2-l) or ^^C0(2-l) transi- 
tions only have upper limits, their line widths are assumed 
to be the same as those of the CO (1-0) and CO (1-0) 
transitions, respectively. When both HCN and HCO''' are 
undetected, we have not attempted to calculate the corre- 
sponding line ratio nor to identify any best model, as the 
ratio of two upper limits is unconstrained. For those cases, 
we do not provide any line width either (see Table [2J . 



3 NON-LTE MODELS 
3.1 Generalities 



As s h own theoretically (e . g. iMeiierink. Spaans fc Israeli 
20071 : ISpaans fc MeiierinkI l2005l : iBavet et all l2008bl. 



2009bl ) and confirmed o bs ervati ona fiy (e.g. Ilsrael fc Baad 



20031: iGao fc SolqmonI |2004|: iKramer et al" 



iGarci'a-B urillo et al.n2006l : iBavet et all l2008al . l2009a 
lAladro et al.l |201G| ). there are several distinct cold gas 
components in external galaxies, each traced by specific 
molecular tracers. For our modeling, we thus consider a 
simple approach with two gas components: one traced 
by the CO lines (a rather extended and cool phase) and 
another traced by the HCN and HCO"*" lines (expected to 
be warmer and denser). This approach is most appropriate 
given the significant differences between the critical line 
densities of these species, i.e. more than two orders of 
magnitude between HCO+(1-0) and ^^CO(l-O), and three 
orders of magnitude between HCN(l-O) and ^^CO(l-O). 
The gas traced by the low-J CO lines and the gas traced by 
HCO^ and HCN are thus not compatible; mixing these two 
components would not allow us to increase our knowledge 
of the interstellar medium (ISM) properties of ETGs as 
much as the two-component approach adopted. We do 



* See |http://www.iram. fr/IRAMFR/GILDAS/docl /html/class- 
html/class, html 



^ More details are avaible in the GILDAS package for both meth- 
ods. 
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Figure 1. New JCMT CO (3-2) line observations of 8 galaxies from our sample. As the line profiles are not always Gaussian, velocity- 
integrated line intensities were extracted by simply summing the intensities over the range of velocities (indicated by the dashed lines) 
with detected 1^00(1-0) emission (see Paper IV). Although IC1024 is clearly detected in ^^CO(3-2), both the shape and scaling of its 
line profile are highly uncertain due to a pointing error. We thus present the IC1024 data here for completeness only. 



note that HCN(l-O) and HCO+(1-0) do not have exactly 
the same critical line density (~ 1 order of magnitude 
difference), and thus they do not strictly speaking trace the 
same gas component. However, in the rest of the current 
paper we will consider that they do, in the sense that they 
both qualitatively trace a gas component denser than that 
observed with the low-J CO lines. 



In the following, we refrain to estimate gas masses and 
gas fractions, as there are large uncertainties in the values 
of the appropriate conversion factors for ETGs (CO-to-H2 
and HCN-to-H2). In fact, there is currently no estimate of 
the HCO^-to-H[2 conversion factor at all in the literature, 
even for late-type galaxies, so it is clearly not possible to 
estimate the gas mass and gas fraction associated with this 
species. Nevertheless, qualitative empirical estimates of gas 
fractions based on the line ratios of ^^CO, ^"^CO, HON and 
HCO^ were discussed in Paper XI. 



3.2 Model description 

To determine the average physical properties of the CO gas 
component on the one hand, and those of the HCN and 
HCO^ gas component on the ot her hand, we use RAD EX , 
the non-LTE model developed bv lvan der Tak et al.l |20o3). 
This model is a zero-dimensional, one molecular species 
code. It solves the radiative transfer equations by assum- 
ing a homogeneous medium and retaining the assumption 
of local excitation. To partially decoupl e the level popula- 
tion a nd radiative transfer equations, the Ivan der Tak et al.l 
i2mi) model uses the escape probability method (one vari- 
ation of which is the LVG method) assuming a uniform (i.e. 
non-expandin g) sphere. T he escape probability equation is 
then given by lOsterbrockl |l974). The model also assumes 
that the optical depth is independent of velocity, appropriate 
for modelling velocity-integrated line intensities. 

Although the assumption of uniform physical conditions 
for each gas component over the scale of a galaxy is crude 
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(e.g. the ISM of nearby galaxies is known to be clumpy), 
these models do provide first order constraints on the aver- 
age physical conditions, and they do represent a useful first 
step in predicting the intensities of the higher order transi- 
tions of CO, HCN and HCO^, on which there is currently 
no constraint in ETGs. In any case, we currently lack infor- 
mation on the very small-scale (< 70 pc) gas distribution of 
ETGs, preventing us from studying their gas properties in 
greater detail. 

Molecular coUisional rates are required to derive the 
molecular gas physical conditions using this model, in addi- 
tion to the spectroscopic and dipole moment dat4j needed 
when using more basic (e.g. LTE) models. We used col- 
lisional rates for CO, HCN and 1100"*" from the Lei- 
den Atomi c and Molecular D a tabas e (LAMDA jj] as recom- 
mended in Ivan der Tak et all (j2007l ). 

There are five main parameters in the LVG model we 
use: the column density N{X) of the species studied X (here 
either ^^CO, "CO, HCN or HCO+) in cm'^ , the line width 
Aw of each transition (see Table [2]) in km s^^, the molec- 
ular hydrogen volume density n(H2) in cm~"^, the kinetic 
temperature Tk in K and the abundance ratios (e.g. iso- 
topic ratios such as ^^C/^C and the HCN/HCO+ abun- 
dance ratio). In our work, we leave the three parameters 
Tk, ?^(H2) and N{X) free (and solve for them by fitting the 
data), whereas we fix the other parameters (Av and isotopic 
abundance ratios) . We note that we do not implement more 
specific characteristics of ETGs in our LVG models, such as 
high metallicity or enhanced a-element abundances. These 
particularities and their detailed influence on the molec- 
ular chemis t ry of ETGs have however been considered in 
iBavet et al.l (|2012l ). 



3.3 Choice of parameters 



As discussed above fSection l2.3|l . the fixed parameter Av is 
taken as the FWHM of the fit for Gaussian profiles or the line 
width of the fit for double-peaked profiles, and is measured 
for each molecular line in each source (see Table El . It varies 
from 58 to 425 km s~^, and its value directly impacts on the 
predicted velocity-integrated line intensity calculated by the 
LVG model. Indeed, the predicted velocity-integrated inten- 
sities are obtained by multiplying the calculated radiation 
temperature of the spectral line by 1.06 An, where the fac- 
tor 1.06 is due to the assumption of a Gaussian line profile. 
Although not all line profiles are Gaussian, the models are 
constrained only by ratios of velocity-integrated line intensi- 
ties, so this factor and the assumed shape of the line profiles 
have no infiuence on the model fits; only variations in the 
line profiles across the different lines involved in the ratios 
matter. This is no longer true for the calculation of the pre- 
dicted spectral line energy distributions (SLEDs), however, 
and the SLEDs of galaxies with non-Gaussian line profiles 
are therefore more uncertain (see Section [4. 3p . 

Currently, there is little knowledge of the true isotopic 
abundance ratio of carbon in ETGs and of its variation from 
source to source. We have thus assumed that it is constant 



^ Prom The Cologne Database for Molecular Spectroscopy; see 
|http: / /w ww . astro. uni-koeln.de/cdms/~ 
' http: / / www. strw. leidenuniv. nl / ~moldata / 



across our sample. As exp l ained for instance bv lMartm et al.l 
(|2010l ) and lRoberts" et al.l (|201ll ) , velocity- integrated line in- 
tensity ratios can be misleading and are not a good proxy for 
abundance ratios, as they are affected by geometry and opti- 
cal depth effects. The velocity-integrated line intensity ratios 
of ^^CO/"CO provided in Paper XI (see their Figure 3) 
can not therefore be used to derive an accurate value of 
the intrinsic ^^C /^^C abundance r atio. S imilarly for HCN(1- 
0)/HCO+(1-0). Iwilson fc RoodI l|l994h reported values of 
the ^^C/^^C abundance ratio in the centre of the Milky Way, 
the 4 kpc molecular ring, the local ISM and the Solar system 
of ~ 20, 50, 70 and 90, respectively. For local st arburst galax:- 
ies, it has been recently estimated to be > 46 (iMartin et al.l 
I2OIOI ). Based on optically thin lines, iRoberts et al.l (|201ll ) 
reported a HCN/HCO^ abundance ratio between 0.5 and 
1.03 for local and extragalactic active star-forming regions. 
For the Galactic centre, only optically thick lines of HCN 
and HCO"'' are present in the literature (and are thus mis- 
leading). Lacking better information for either case, we have 
therefore adopted a fixed value of 70 for the ^^C/^''C abun- 
dance ratio and a value of 1 for HCN/HCO^. Only knowl- 
edge of the small-scale distribution of the ^^CO, "CO, HCN 
and HCO"^ emission for various line transitions would allow 
us to properly estimate the true abundance ratios. 

We thus built four grids of LVG models to be combined, 
one for each of ^^CO, "CO, HCN and HCO+. Each grid 
covers kinetic temperatures of 10 to 250 K in steps of 10 K, 
gas volume densities of 10^ to 10^ cm"'' in steps of 0.5 dex, 
column densities of 10^" to lO'^" cm"'^ in steps of 0.5 dex 
and line widths of 60 to 430 km s^^ in steps of 10 km s~^. 
Each grid thus contains over 150, 000 models, leading to a 
total of over 600, 000 models for the species investigated in 
the present study. 

The model outputs are numerous, and among them is 
the quantity we are primarily interested in, i.e. the predicted 
velocity-integrated line intensity (in K km s"'^). We ran each 
model to include the first 15 energy levels and therefore the 
first 15 transitions of each species, allowing us to predict 
the most interesting observational lines for Herschel Obser- 
vatory and the Atacama Large Millimeter/sub-millimeter 
Array (ALMA) (see Section g^J. 

For each galaxy, we have thus compared the predicted 
velocity-integrated line intensity ratios (computed directly 
from the model outputs) with the observed ones. More 
specifically, we have worked with the independent ^^CO(l- 
0)/"CO(1-0), '^CO(2-l)/"CO(2-l), i2CO(l-0)/'2CO(2-l) 
and ^^00(1-0)7^^00(3-2) (whenever the ^^CO(3-2) line was 
available; see Table [TJ line ratios on the one hand, and the 
HCN(l-O) /HCO+(1-0) line ratio on the other hand. For each 
ratio, we have corrected for beam dilution as described in 
Paper XI. The ^^00(3-2) line emission is only used in one 
ratio, not two, because we wanted to keep independent con- 
straints for reliable calculations (see Section [4. 2|l . 

To identify the model whose input parameters (Tk, 
n(H2) and N{X)) reproduce best the observed line ratios, 
we have therefore constrained the models with a minimum 
of three observed line ratios for the CO gas component, and 
only one observed line ratio for the HCN and HCO^ gas 
component. The physical properties derived for the CO gas 
component are thus more reliable than those derived for the 
HCN and HCO"'' one, as the latter are under-constrained. 
The results for the HCN and HCO"'' gas component should 
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Table 2. Line widths for our sample of 18 gas-rich ETGs, extracted from the observations published in Paper XI and from the new 
^■^CO(3-2) observations presented here. 



Galaxy D At; (12CO(1-0)) Ai; (12C0{2-1)) At; (i2cO{3-2)) At; (i='CO(l-0)) At; (13C0(2-1)) At; (HCN(l-O)) At; (HCO+ (1-0)) 

(Mpc) (km s-i) (km s-^) (km s'^) (km s'^) (km s-^) (km s-^) (km s-^) 



IC0676 


24.6 


97.3 ± 17.4 


87, 


.1 ± 10.7 


84.1 ± 3.1 


150.2 ±39.1 


77.7 ± 16.1 


101.2 ±28.1 


116.8 ± 26.7 


IC1024 


24.2 


160.9 ±19.3 


146, 


.9 ±22.4 


123.1 ± 20.0 


191.0 ±30.6 


152.6 ± 19.8 


- 


- 


NGC1222 


33.3 


158.7 ± 10.4 


169, 


.8 ± 15.7 




170.8 ±31.8 


133.7 ± 17.7 




128.7 ±30.0 


NGC1266 


29.9 


128.1 ±22.3 


152, 


.2 ± 28.9 




95.1 ± 19.1 


104.0 ± 20.5 


157.3 ±48.3 


179.6 ± 13.7 


NGC2764 


39.6 


249.5 ± 6.6 


226, 


.0± 7.0 




250.9 ±22.8 


246.3 ± 25.0 


220.3 ±62.4 


199.2 ±34.6 


NGC3032 


21.4 


128.8 ± 14.5 


102, 


.5± 3.3 




80.8 ± 16.5 


103.4 ± 11.2 


162.5 ±45.5 




NGC3607 


22.2 


268.6 ±40.5 


323 


.1± 4.7 




425.9 ± 85.2 


386.7 ±68.6 


385.4 ±71.4 


387.3±77.8 


NGC3665 


33.1 


318.1 ±64.7 


319 


.4 ±51.7 


317.8 ±54.7 


327.5 ±65.5 


312.6 ± 76.8 


318.3 ±27.8 




NGC4150 


13.4 


159.2 ± 13.8 


189, 


.1 ± 17.6 




149.0 ± 29.8 


158.2 ±30.4 






NGC4459 


16.1 


194.6 ± 35.7 


192, 


.7 ±28.9 




196.3 ±33.3 


190.1 ±39.3 


193.3 ±42.5 




NGC4526 


16.4 


348.9 ± 34.9 


341, 


.1±68.0 


336.7 ± 53.8 


345.6 ±55.3 


345.2 ±72.5 


346.7±73.4 




NGC4694 


16.5 


80.2 ± 5.2 


70, 


.9 ± 14.2 




57.9 ± 9.3 


100.8 ±30.3 






NGC4710 


16.5 


229.2 ±45.3 


209, 


.5 ±42.0 




292.7 ±52.8 


261.8 ±52.1 


243.0 ±59.5 


272.5 ±43.4 


NGC5866 


14.9 


294.8 ± 50.5 


242, 


.5 ±48.6 


240.7 ±43.3 


374.3 ± 71.1 


377.9 ±74.1 


338.6 ±71.9 


360.5 ±33.6 


NGC6014 


35.8 


179.7 ±28.8 


122, 


.1 ± 24.6 


107.2 ± 14.0 


162.1 ±25.9 


140.8 ±29.5 


93.3 ± 17.9 




NGC7465 


29.3 


129.8 ±24.8 


80, 


.1 ± 12.7 




140.8 ± 19.7 


130.9 ±27.6 


150.5 ±42.6 


124.7 ± 12.9 


PGC058114 


23.8 


155.8 ± 29.6 


160, 


.5 ±32.1 


149.4 ±22.4 


174.9 ± 27.9 


154.7 ±32.5 




170.9±54.4 


UGC09519 


27.6 


179.5 ±32.1 


155, 


.0±31.5 


128.0 ± 15.7 


153.1 ±26.0 


199.5 ±46.9 




141.2±34.2 



Notes: The symbol '-' implies that no line width could be measured, either because there is no observation (e.g. ^-^00(3-2) line) or 
because signal-to-noise ratio of the line is too low. However, when required, the line width adopted is that from the 'closest' species (see 

Section [2} . 



therefore be considered purely as indicative, and we present 
these results (table and figures) in Appendix|A]only. We will 
focus exclusively on the CO gas component for the remain- 
der of the current paper. 

Finally, we note that the observed line ratios in- 
volving isotopologues (^2CO(1-0)/^^CO(1-0) and ^^00(2- 
l)/^''CO(2-l)) have been reproduced by combining a ^^CO 
model having a column density A'^(^^CO) with a ^^CO model 
having a column density iV("CO) = (^^C/^2C)7V(^^CO). 
Similarly, an observed ratio line involving a ^^CO(l-O) line 
with a FWHM of 300 km s"^ and a ^^C0(2-l) fine with a 
FWHM of 290 km s~^ has been reproduced by combining a 
^^CO model having a velocity width of 300 km s~^ for the 
(1-0) transition with a ^^CO model having a velocity width 
of 290 km s~^ for the (2-1) transition. 



4 RESULTS 

We explore below two complementary methods to identify 
the 'best' model for each galaxy studied: the best-fit model 
based on the method (see Section[4T]) and the most likely 
model based on the likelihood (see Section [4. 2p . From these 
two methods, we calculate for the first time the correspond- 
ing predicted spectral line energy distributions of CO in local 
ETGs (see Section gS]) . 

4.1 Best-fit model identification 

For each source, we identify a best model using the min- 
imisation method for each of the two gas components inves- 
tigated, i.e. the gas traced by CO on the one hand and the 
gas traced by HCN and HCO^ on the other. The physical 
properties (Tk, n.(H2) and N{X)) corresponding to these 
best fits are listed respectively in Table [3] for the CO gas 



component and in Table IATI for the HCN and HCO"*" com- 
ponent. This is the first time that the physical properties of 
the molecular gas are constrained in ETGs in a systematic 
and homogeneous way, and for a relatively large number of 
sources (18). These properties are of prime interest for fur- 
ther individual and statistical studies of the gas in ETGs 
(physical conditions, origin, Hl-to-H2 ratio, etc). When cal- 
culating the x^5 the errors on the observed line ratios have 
been properly calculated and are fully taken into account. 
As the errors on line ratios are asymmetric (see Paper XI), 
we have considered the geometric average of the two values 
in our x^ calculations. 

The kinetic temperatures obtained for the best-fit mod- 
els are either low (Tk = 10-20 K for IC1024, NGC1222, 
NGC1266, NGC2764, NGC3032, NGC3665, NGC4150, 
NGC4459, NGC4526, NGC4694, NGC4710 and NGC5866), 
intermediate (30-70 K for IC0676, NGC3607, PGC058114 
and UGC09519) or high (120 K for NGC6014 and 150 K 
for NGC7465). Although these kinematic temperatures may 
seem high in view of the commonly assumed low star- 
formation activity of ETGs, the fact is that the absolute 
star-formation rates (SFRs) of ETGs are not always low (es- 
pecially in the molecular gas-rich objects considered here). 
It is the specific SFRs (SSFRs; i.e. the SFRs per unit stellar 
mass ) that are generally lower than those of spirals (see, 
e.g., IShapiro et al] l20ld for the SAURON ETGs). ETGs 
also have additional heating, excitation and ionisation mech- 
anisms that are unlikely to be important in other galaxy 
types. For example, old (and therefore not massive) but hot 
stars will heat the dust grains and, through the photoelec- 
tric effect, the gas, impacting both the molecular and ionised 
gas emission lines. This peculiar population als o gives rise 
to the ultraviolet (UV)-upturn phenomenon (e.g. lO'Connelll 
ll999l : lBureau et al.llioil] : Ijeong et al.ll2012l ). increasing the 
UV flux (and thus the heating) surrounding the gas indepen- 
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dently of star formation. Supernova rates leading to cosmic- 
ray fluxes may also be different than those of spirals. Finally, 
the high stellar metallicity and g-ele ment enha nc emen t may 
also play a role in the gas chemistry (jBavet et al.ii2012l ') . Any 
or all of these mechanisms may thus partly be responsible 
for the high gas temperatures measured in some galaxies, 
although it is impossible at this point to determine which 
has the greatest influence. 

The best-fit H2 volume densities do not vary by 
more than a factor of 10 between our sources (n(H2) = 
10^-* cm-3), except for NGC1222 and NGC5866 that have 
n(H2)= 10®-^ and 10^ cm"^, respectively. Similarly, the CO 
column densities derived are pretty constant amongst all 
the galaxies studied, with N{CO)^ ]^q18.5-20 ^^-2 rj-j^ggg 
column densities may appear high, especially once con- 
verted to N(H2). Indeed, assuming a standard Milky Way 
CO-to-H2 abundance ratio of 10~*, one obtains N(H2) = 
2^022.5-24 j,jjj-2^ coresponding to an average surface mass 
density of 10^'^~* Mq pc~^. For comparison, giant molec- 
ular clouds (GMCs) in the Milky Way have ty pical surface 
densit ies of 10^ Mq pc~^. However, as shown bv lBavet et al.l 
who studied in details the chemistry of ETG-like 
environments, the standard Milky Way CO/H2 fractional 
abundance is not appropriate for ETGs, and a value highe r 
than 10"'' should be used (see Fig. 2 in lBavet et"alll2012h . 
For most of the sources studied here, a value of 10~^'^ is 
more appropriate, leading to surface densities closer to the 
Milky Way value. In addition, high CO column densities may 
not be totally unexpected, as sugg ested by our CO interfer- 
ometric maps lAlatalo et al.ll2012l . hereafter Paper XVIII), 
where high surface densities are observed despite low SS- 
FRs (as high as ^ 10* Mq pc"^ in NGC1266). These CO 
column densities may also be considered as lower limits, 
since the authors conservatively assumed a standard CO- 
to-H2 conversion factor. This apparent low star-formation 
efficiency (SFE; i.e. the star formation rate per unit dense 
gas mass) may partly be due to the different dynamics of 
ETGs compared to that o f spirals (morphological quench- 
mg; see iMartig et al1l2009l V Nevertheless, most ETGs show 
typical surface densities of « 10^ Mq pc~^^, consistent with 
Milky Way values. 

Table [3] reveals that the galaxies detected in ^^C0(3- 
2) are not systematically those showing the highest kinetic 
temperatures and gas densities as determined by the best- 
fit method (see, e.g., NGC3665 and NGC4526). As it takes 
warmer and denser gas to effectively populate the (3-2) level, 
this may seem surprising at first, but in fact our best-fit 
results only show that the CO (3-2) line in these galaxies 
is not the main driver of our best-fit model identification. 
In other words, there is basically nothing special about the 
objects detected in ^^^00(3-2). In fact, we did not select for 
observations objects predicted to have the highest absolute 
^^CO(3-2) fluxes, but rather the highest apparent ^^CO(3-2) 
fluxes (e.g. because they are nearby or have high ^^CO(l-O) 
fluxes). So the other objects are not so much non-detections 
as they are non- observations (we just do not know what their 
^^CO(3-2) fluxes are, and they may well be high). A larger 
sample of CO (3-2) observations in ETGs is thus needed 
for further comparison. 

Understanding the exact driver of the best-flt model 
identiflcation requires a detailed analysis of the LVG code 
outputs and of the minimisation, that are outside the 



scope of this paper. However, to evaluate the accuracy of 
our best-flt CO models and calculate uncertainties on the 
derived best-flt model parameters listed in Table O we have 
calculated the A^^ = ~ Xmin contours for each source 
in our 3D parameter space (Tk, n(H2), A'^(CO)). The flrst 
flve contours in Figure [5] show the 1 to 5o[f| confldence level 
contours in three TK-n(H2) planes centred on the best-flt 
A'"(CO) value, where dark-grey dots (red dots in the online 
version) correspond to individual models. The Ax^ contours 
for the HCN and HCO^ gas component are presented in Fig- 
ure |Xl] Additional contours of Ax^ =50, 100, 200 and 300 
are used in these figures only for guidding the eyes further 
and see better how the confidence contours are spreading 
within the model grid. 

With additional (higher J-levels) observations of the 
same species, probing denser and warmer gas closer to the 
forming stars, we should be able to «^ confirm the results 
obtained here and ii) understand better the b est-fit model 
parameter differences across galaxies (see, e.g., iBavet et al.l 
|2004 12OO6I who used the ^^00(6-5) and ^^00(7-6) lines in 
late- type galaxies). 

Despite the uniqueness of the best-fit model for each 
galaxy (light-grey filled circles in Fig. (2] green filled circles 
in the online version), there is a range of models included in 
the Icr confidence level (darkest zones in Fig. [2} . These mod- 
els can all also be considered as good models, and they could 
be used to define formal uncertainties, although these would 
be artificially large along the three axes simultaneously and 
would not convey the complex shape of the Ax^ contours. 
Indeed, as seen in Figure [2] the contours generally have a 
banana-like s hape similar to w hat is often obtained in disc 
galaxies (e.g. iBavet et al.ll2006l ). even if the contours are of- 
ten more ragged (explained by the fact that there are fewer 
contraints in our study). This means that the kinetic tem- 
perature and gas volume density solutions are degenerate. 
Despite this, we note that the low a confidence level contours 
are relatively narrow in the 3D parameter space for IC0676, 
NGC3032, NGC4526, NGC5866 and UGC09519. The uncer- 
tainties on the best-fit parameters (Tk, n(H2), A'^(CO)) are 
thus relatively small and the best-fit models identified are 
more reliable for these sources. 

4.2 Most likely model identification 

To better understand the distribution of good models in our 
3D parameter space, and to define robust uncertainties, we 
have investigated a complementary 'best' model identifica- 
tion method that consists in determining the marginalised 
probability distribution functions (PDFs) of the three in- 
put model parameters (Tk, n(H2) and A ^ (CO)) , following 
the formalism developed in iKavirai et al.| (|2007l ). For each 
parameter value, we have calculated the likelihood, i.e. the 
sum of the exponentials of — Ax^/2, marginalizing over the 
other two free parameters. This technique allows us to iso- 
late the effect of each parameter, and the results are shown 
in Figure [3] 

From the plots shown in Figure |3] we have calculated 
the total area defined by each PDF, terminating the integra- 
tion when necessary at the end of our model grid. This area 

* The <T symbol represents the standard deviation. 
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Figure 2. Best-fit models for the CO gas component. Each panel shows the Ax^ 



contours of the CO gas component as a 



function of the kinetic temperature Tk and H2 volume density n(H2), for three values of the CO column density N{CO) (hence three 
plots) centered on the best-fit value and indicated at the bottom of each plot. The model grid is shown with dark-grey dots (red dots 
in the online version), whereas the best- fit model of each galaxy (x^ • or Ax^ = 0) is labeled with a light-grey filled circle (green filled 
circle in the online version). The Ax^ contours and greyscales are for 1 to 5a confidence levels (i.e. Ax^ = 3.54, 6.25, 7.81, 9.35 and 
11.34 when three line ratios are used and Ax^ = 4.73, 7.78, 9.49, 11.14 and 13.27 when four line ratios are used). Additional contours 
of Ax^ =50, 100, 200 and 300 are used in these figures only for guidding the eyes further and see better how the confidence contours 
are spreading within the model grid. The low a, high confidence level contours containing 'good' models are represented by the darkest 
areas. The galaxy name is indicated in the top-right corner of each panel. 
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Table 3. Best-fit model parameters of the CO gas component. 



Table 4. Mean likelihood parameters and their la uncertainties. 



Galaxy 




Tk 


"(Ha) 


Af(CO) 






(K) 


(cm~^) 


(cm~2) 


lOUU 1 \> 


1 n7 


1 u 


1 o3 


ini9.5 

-LU 




U.OD 




1 0* 




IVr'PI 999 




9n 
zu 


1 0^.5 


10l9 

-LU 




u.oo 


90 
zu 


1 0* 


inl8.5 
-LU 


IN VT V_.' Z ( U'i 


0.11 


10 


igS.S 


IQlS 


NGC3032 


5.32 


10 


103.5 


1019 


NGC3607 


0.23 


50 


10^ 


1020 


NGC3665=» 


6.45 


10 


103.5 


1020 


NGC4150 


0.02 


20 


103.5 


1019 


NGC4459 


0.57 


10 


10* 


1019.5 


NGC4526'' 


1.28 


20 


103.5 


1020 


NGC4694 


2.15 


10 


103.5 


1018.5 


NGC4710 


1.38 


20 


103.5 


1019.5 


NGC5866'' 


41.18 


10 


10^ 


1019.5 


NGC6014'' 


1.45 


120 


10^ 


1020 


NGC7465 


9.95 


150 


104.5 


1020 


PGC058114'' 


0.34 


30 


103.5 


1019 


UGC09519'' 


25.21 


40 


10^ 


1019 



Galaxy 



(K) 



n(H2) 
(cm"-^) 



Af(CO) 
(cm-2) 



Notes: ^: Galaxies with a i2cO(3-2) line flux used in the mod- 
elling work. 



represents the total probability (100%) of each parameter, 
from which we have determined the mean. This defines the 
equipartition (50%-50%) of the PDF and, for well-behaved 
(i.e. single-peaked and symmetric) PDFs, is usually consid- 
ered the most natural way to identify the most likely model. 
The mean is shown in each plot of Figure [3] by a solid grey 
line surrounded by two dashed grey lines, themselves defin- 
ing the Icr uncertainty (68% probability, 34% on each side 
of the mean). The mean and associated la uncertainties are 
listed in Table |4] for each model parameter and galaxy. 

Similarly to the best-fit method, the likelihood method 
yields groups of galaxies with similar ranges of kinetic tem- 
peratures. These categories are however less marked than for 
the best-fit models (see Section BTTl and Table[3}. The H2 vol- 
ume densities and CO column densities are again rather uni- 
form (variations of less than a factor of 10 for most sources, 
excluding NGC1222 and NGC5866), similar to what is ob- 
served with the best-fit method. 

The likelihood method helps us to identify the 'best' 
model in a complementary way to the simple method. 
Comparing Tables [3] and |4] and the positions of the black 
dots (best-fit models) and solid grey lines (most likely mod- 
els) in Figure |3l we see that in the majority of sources 
the two methods yield consistent results within the (admit- 
tedly rather large) uncertainties. Indeed, one notes that the 
highest PDF peaks agree well with the best-fit model para- 
maters. In a few cases (NGC1266, NGC6014 and NGC7465 
are the best examples), however, the most likely models do 
not exactly correspond to the best-fit models. The apparent 
discrepancy between the best-fit and most likely models is 
then generally due to the double-peaked (or triple-peaked 
for NGC7465) shape of the kinetic temperature PDF (see 
the left panels in Fig. [Sjl . These two (three) peaks, when of 
almost equal amplitudes, imply that there are two (three) al- 
most equally probable solutions. The derived mean for these 
distributions, located at a roughly equal distance from the 
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Notes: Sources with a marginalised probability distribution 
function with two or more peaks, where the mean parameter listed 
in the table is misleading, i.e. it is not the most likely parameter 
(see text in Section l4.2l l. 

two (three) peaks, is therefore misleading and does not gen- 
erally correspond to the most likely model. 



4.3 ETG CO spectral line energy distributions 

Calculating models up to the 15*'' level population for 
the CO molecule allows us to predict velocity-integrated 
line intensities up to i'^CO(15-14). Converting the velocity- 
integrated line intensities (in K km s~^) into int ensities (in 
W m"^ sr"i) using Eq. 4 in iBavet etall (|2004l '). and plot- 
ting both quantities as a function of the higher rotational 
level of the transition considered (hereafter Juppcr), yields 
the predicted CO SLED of our ETGs (see Fig. [J). 

In Figure m we show the best-fit model SLED as a solid 
black line and use two solid grey lines to delineate the range 
of SLEDs associated with a la confidence level on the model 
parameters. In other words, all the models included within 
the first Ax^ contour in Figure [5] have their SLEDs between 
the two solid grey lines. The SLED corresponding to the 
most likely model is also shown as a dotted black line. Unsur- 
prisingly, the largest la SLED ranges occur for sources with 
a large kinetic temperature uncertainty (i.e. when the 3D 
darkest region in Fig. [5] has a large extent along the y axis). 
For example, within the la confidence level of its model 
parameters (see Fig. [Sjl , the predicted CO SLED peak posi- 
tion (hereafter J^ax) of NGC6014 varies between Juppcr — 5 
and 7, covering three increments in Juppcr, while its kinetic 
temperature varies by a factor of 6. A similar behaviour 
is observed for NGC1266, NGC4694, NGC4710, NGC7465 
and PGC058114, with a predicted SLED turnover position 
increment of 2 or 3 and a kinetic temperature increase by 
a factor of 3 to 6 within la. When Tkin does not vary by 
more than a factor of 3, Jmax only varies by one increment 
in Juppcr or less (e.g. IC0676 and IC1024). This provides a 
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Figure 3. Marginalised probability distribution functions of our model parameters. For each galaxy (indicated in the top-right corner of 
the left-hand plots), the PDF of each model parameter is shown, marginalised over the two others. The best-fit (x^i^) model parameters 
are indicated by solid black circles on the x axes. The mean likelihood values and Icr uncertainties are represented by solid and dashed 
grey lines, respectively. 



clear indication of how the peak of a galaxy SLED varies 
with kinematic temperature, and appears consistent with 
other current theoretical studies (e.g. Wu et al., in prep.). 
Of course, the H2 volume density also has an important im- 
pact on the shape of the SLEDs, but with our present galaxy 
sample the corresponding effect is less clear (see Fig.|4]). It is 
however generally thought that when the gas volume density 
increases, the SLED widens (Wu et al. in prep.). 



To our knowledge, ours are the fir st CO SLEDs ever 



produced for normal ETGs (but see Eckart et al. 1990l. 
'Wild. Eck art fc Wiklindlll997l and iMuller fc Dinh-V-Trund 
2009 for modelling results on Centaurus A). They are thus 
of prime importance to better understand the star forma- 
tion activity of these systems. Indeed, Jmax can be used as 
a tracer of a galeixy's star formation activity, as it is this ac- 
tivity that directly drives the gas excitation in most sources. 
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including those studied here. As seen in IWeifi et all (|2007l ) 
for late-type galaxies, the position of the SLED turnover 
appears linked to the star formation rate. For example, the 
galaxy APM0879 (SFR> 100 Mq yr"^) has the highest 
observed SLED turnover position, whereas the Milky Way 
(SFRss 1 Mq yr~^) has the lowest. It has also been noted 
that the SLEDs of nearby (< 10 Mpc) starburst galax;ies 
peak at higher transitions (around (7-6)) on average than 
those of nearby normal spirals (such as the Milky Way, 
whose SLED turnover is located at around the (4-3) transi- 



tion; see, e.g., iBavet et al.ll2004 l200d ). Nevertheless, there 
is currently no study showing a direct link between a SLED 
turnover position and the SFR. What drives the position 
of the turnover is thought to be not the kinematic tem- 
perature alone, but rather a combination of the kinematic 
temperature and gas volume density. The gas column den- 
sity is expected to impact mostly on the intensity of the 
SLED turnover, not its position. The ^^C/^"^C fractional 
abundance (for CO), line width A« and opacity can also in- 
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fluence the SLED turnover position, but they are secondary 
order effects (again, see Wu et al. in prep.). 

We present in Table[4]a summary of the predicted SLED 
turnover positions (Jmax) for our galaxies, derived from the 
best-fit models (best-fit and Icr uncertainties). We also show 
in Figure [5] (top) those SLED turnover positions as a func- 
tion of each of the best-fit model parameters (Tk, n{]l2) and 
A'^(CO)). There are weak trends of Jmax with both Tk and 
A''(CO). Those trends however need to be confirmed with 
a more reliable (i.e. larger) statistical sample, as they are 



mostly driven by the modelling results of a single galaxy 
(NGC7465, with the highest T^in and Jmax). In addition, 
with only the (1-0), (2-1) and occasionally the (3-2) tran- 
sition, it is often difficult to reliably predict the existence 
and/or position of a SLED peak, and therefore to confi- 
dently infer the level of star formation activity. The results 
presented here, although suggestive and exciting since avail- 
able for the first time, are thus based on limited data and 
a simple theoretical model with many limitations and as- 
sumptions, and should therefore be considered accordingly. 



16 Bayet et al. 





50 100 150 200 250 



^(GC6014^ 



lO'" 



50 100 150 200 250 



O 

o 

2 



3 




I I I I I I I 

NaC7465 



1 .1 1 



50 100 150 200 250 
Tk (K) 



10 



10 



5 
4 
3 
2 
1 






1 ' 












6 - 


1 






4 - 


: 0.5 














2 - 





1 




n - 



□ 10 






1 


6 r 


1 1 1 1 1 

CJ CO 








4 - 


: 4 


;1 
















2 - 


: 2 










= 








- 


- 11 


T' " T " t 



10 



18 



10 



20 





n(H2) (cm 

Figure 3. (Continued). 



N(CO) (cm ^) 



Recently, IShapiro et all (|20ld ) and Falcon-Barroso et 
al. (in prep.) derived preliminary star formation rates for 8 
of the molecular gas-rich ETGs in our study (NGC1266, 
NGC3032, NGC3607, NGC4150, NGC4459, NGC4526, 
NGC4694 and NGC5866), with SFRs ranging from 0.06 
to 0.59 Mq yr~^ (see Table [5]). These estimates are based 
on non-stellar 8 ^m emission. Obtaining SFR estimates for 
ETGs is however fraught with danger, as most usual SFR 
indicators are not straightforwardly applicable to ETGs. At 
the very least, specific properties of ETGs (e.g. UV-upturn 



population. X-ray halos, etc) may significantly alter the cal- 
ib ration factors used in those methods (see the discussion 
in lCrocker et aLlbOllh . Figure [5] (bottom) shows the SLED 
turnover position as a function of the 8 /.tm-derived SFR for 
those 8 galaxies. No clear trend is observed. Plots showing 
these SFR estimates as a function of the best-fit model pa- 
rameters (Tk, n(H2) and A''(CO)) show no clear trend either 
(Fig. [S| . The sample size is however rather small, prevent- 
ing us from drawing reliable conclusions. Whether or not 
ETG SLED turnover positions follow the same trend with 
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SFR as those of late- t ype galaxies (e.g. IWeifi et all l2007l : 
iBavet et al.]|2004 [ 20061 ) thus remains to be confirmed. Be- 
cause of the complicating issues with ETG SFRs mentioned 
above, a detailed study of the SFRs of our sample galaxie 
is beyond the scope of this paper, and we refer the reader 
to lShapiro et al.l (|2010l ) and Falcon-Barroso et al. (in prep.) 
for further discussions. 

Despite their limitations, our models are nevertheless 
very useful tools, making specific predictions and paving 
the way for future observations. These predictions can in 
turn inform on the models themselves. For instance, the two 
methods of 'best' model identification yield different pre- 
dicted turnover positions in some sources. If the ^^C0(6- 
5) line was detected in IC1024 and NGC1266, for exam- 
ple, with fiux in both cases below 10"^" W m~^ sr~^, a 
flat SLED seems an unlikely possibility for those two galax- 
ies. In such case, there should rather be a SLED turnover 
and it should be located between J=2-l and J=6-5. On the 
contrary, if the ^'^ CO (6-5) line was detected in these two 
galaxies but with a flux of about 10~^ W m~^ sr~^, this 
will favour a flat SLED scenario thereby indicating that 
more sophisticated models are required to infer the prop- 
erties of the molecular gas (e.g. photon-dominated or cos- 
mic ray-dominated region codes such as those developed by 
Le Petit ct al. 2006:'Bavet et al.ll2009bl . l201ll :|PapadopouIos' 
2010; Mciicrink ct al. 2011,). Finally, the great advantage of 
detecting the ^'^ CO (6-5) line, in most sources displayed in 
Figures U would be to single out the most accurate method 
of best model identification in ETGs (either or maximum 
likelihood) . 

Further limits to the accuracy of the predicted SLEDs 
arise from the underlying assumptions, that have to be 
carefully considered. For example, the shapes of the line 



profiles are neglected when working with line ratios (see 
Section 13. 3p . However, this is not true anymore when 
modelling velocity-integrated line intensities or intensities, 
and unfortunately RADEX assumes a Gaussian line profile 
(j yan dor Tak et al. 2007). We know that for galaxies such 
as NGC3032, NGC3607, NGC3665, NGC4459, NGC4526, 
NGC5866 and NGC6014, a double-peaked line profile is 
more realistic, so for these sources the SLEDs shown in Fig- 
ure |3] are more uncertain. 



5 DISCUSSION 

5.1 Correlations with global parameters 

In Paper XI, ratios of molecular lines discussed 
in this paper were shown to correlate with some 
global physical parameters of the sample galaxies 
(e.g. absolute A's-band magnitude, Hl-to-H2 mass ra- 
tio, dust morphology, etc). It is thus natural to ask if 
the molecular gas physical conditions derived here 
show similar correlations. We therefore investigated if 
the rough groups of sources identified in Sections 14.11 and 
14.21 (see Tables |3] and |4]) , based on common molecular gas 
physical conditions, are related to other characteristics of 
the galaxies derived from various other ATLAS'^'^ datasets. 
We investigated potential correlations of our derived gas pa- 
rameters with most available quantities: distance, absolute 
T^s-band magnitude, kinematic misalignment angle of the 
molecular gas with respect to the stars, inclination, HI mass, 
H2 mass, 60 to 100 fim flux ratio (/eo/Zioo), H/? emission 
fine equivalent width, H/3, Mg6, Fe5015 and Fe5270s absorp- 
tion line strength indices, stellar age, metallicity ([Fe/H]) 
and a-element (over-)abundance ([a/Fe]), and effective ra- 
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Figure 4. Predicted ^^CO SLEDs. For each galaxy (indicated in the top-right corner of each plot), the SLED is expressed in K km s~^ 
(left) and in W sr~^ (right). The solid and dotted black lines show the predictions from the best-fit and most likely model, 

respectively. The solid grey lines define the range of possible SLEDs associated with the Icr confidence level on the best-fit model 
parameters (see Section 14.311 . Crosses indicate the model predictions. The ^■^CO observations are represented by solid black circles with 
error bars. 
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Figure 4. (Continued). 
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Figure 4. (Continued). 
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Figure 4. (Continued). 
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dius (Rc). Most of these p arameters were obtained fro m the 
ATLAS^° survey papers l|McDermid et al.l l20 06b'a'. 2007; 
Paper I; Paper II; Paper III; Paper IV: iDavis et al..,2011b„ 
hereafter Paper V; Paper VII; Paper X; Paper XI) or from 
Sarzi et al. (priv. comm.). No correlation has been identi- 
fied with any of the three model parameters (Tk, n(H2) and 
A'^(CO)) reproducing best the observations. 

A possible explanation for the absence of any correla- 
tion is that, although the low-J CO transitions are good 
tracers of the total molecular gas reservoir, only higher-J 
CO transitions closer to forming stars may correlate well 
with, e.g., stellar properties. Or perhaps the objects stud- 
ied here simply do not span a large enough range in stellar 
mass. Only future observations of higher-J CO transitions 
(e.g. ^^CO(6-5)) will confirm such a hypothesis. If a sub- 
stantial fraction of the molecular gas in ETGs has an exter- 
nal origin, as argued in Paper X, then correlations with the 
dominant old stellar population may in fact not be expected. 
Hence, the absence of correlations may simply support the 



case for an external gas origin. However, given the ex- 
istence of correlations with the observed line ratios 
(Paper XI), the absence of correlations with the de- 
rived physical conditions is puzzling. It may there- 
fore also suggest that the models considered here 
are too primitive and do not do justice to the likely 
complexity of the ISM in the sample galaxies. Multi- 
component models mixing photons dominated, X- 
rays dominated and cosmic-rays dominated regions 
will be investigating soon. 



5.2 Literature comparison 

To further understand the star-formation activity of our 
molecular gas-rich ETGs, we compared their SLEDs and 
best-fit model parameters to those obtained for well-known 
nearby galaxies (< 10 Mpc), including starbursts (Maffei2, a 
young starburst; NGC253, an intermediate-stage starburst; 
and M82, a post-starburst), a cosmic ray-dominated galaxy 
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Table 5. Turnover positions of the best-fit model SLEDs. 



Galaxy 


J max 


SFR 
(M0 yr-i) 




5-0 










MP Pi 999 


fi+1 




NPP1 9fifi 


4+2 






o+l 






o+l 
•^-0 




NGC3607 


51? 


0.22 


NGC3665 


3tJ 




NGC4150 


4t? 


0.06 


NGC4459 


4t? 


0.17 


NGC4526 


4+1 


0.37 


NGC4694 


q + 2 


0.17 


NGC4710 






NGC5866 


4+0 


0.12 


NGC6014 






NGC7465 






PGC058114 






UGC09519 


3lJ 





Notes: Jmax uncertainties are defined as the maximum Jmax 
ranges corresponding to a Icr u ncertainty on the best -fit model pa- 
rameters. SFR values are from lstiapiro et al. I ll2010h and Falcon- 
Barroso et al. (in prep.). 

(NGC6946), AGN-dominated galaxies (Markarian231 and 
Arp220) and normal spiral galaxies like the Milky Way. This 
small sample of nearby galaxies is not exhaustive, but it 
rather aims to be representative of the various types of star- 
formation activity present locally. We can thus attempt to 
assess whether the processes at work in these well-known 
'archetypes' of each star-forming class are the same (or give 
rise to the same gas excitation conditions) as those in ETGs. 
We know that the formation and evolution of ETGs and 
late- type galaxies are very different, but are these differ- 
ences also reflected in the star-forming gas (also confined 
to a disc in the majority of our sample galaxies; e.g. Pa- 
per XIV; Paper XVIII)? Only a rough comparison is pos- 
sible here because of the quality and incompleteness of the 
datasets available in the literature, but it is a useful first 
step. 

To reproduce the low-J CO line emission of the 
young starburst galaxy Maffei2, a kinetic temperature 
Tk > 50 K, an H2 volume density n(H2)~ 10^ '^ cm"^ 
and a CO c olumn density A^(CO)« 5-7 x 10^'^ cm"^ 
are required |W cliachc w^ Casoli fc Combeslll988l : iHurt et al] 
I1993I : iDumke et al.i i200lll . A prediction of the CO SLED 
turnover position is not available in the literature. 
iBavet et al.l l|2004l ) showed that the CO SLED of NGC253, 
at an intermediate starburst evolution stage, peaks at 
CO (6-5), whereas their best-fit model parameters are 
Tk = 70-150 K, n(H2)> 1 x 10* cm"^ and N{CO)= 
1.5 ±0.5 X 10^^ cm~^. Recent Herschel Observatory results 
on M82, a post-starburst galaxy, show that the CO SLED 
peaks at ^^CO(7-6) and that Tk = 350-825 K, n(H2) = 
10^"*-^ cm"^ and N{CO)= ipis-s-i g-s cm~^ are require d 
to reproduce the observed CO data (|Panuzzo et al.l [20101 ). 
We note however that observations up to the ^^CO(13-12) 
transition have been obtained, and that several gas com- 



ponents (whose parameters range as above) are needed to 
reproduce the full SLED. The authors do not provide the gas 
parameters for each component, so only the lower values of 
the ranges me ntioned here are ap propriate for comparison 
with our study. iBavet et al.l (|2006t ) showed that the CO gas 
component in the cosmic ray-dominated galaxy NGC6946 
requires Tk = 130 K, n(H2)= 1.5 x 10^ cm"^ and A^(CO) = 
4.4 X 10^" cm"^ to be excited, and that its CO SLED 
peaks at CO (5-4). The AGN-domina ted galaxy Markar- 
ian23 1 shows an unusual CO SLED (|van der Werf et al.l 
l2010h . becoming flat after the CO (6-5) transition. Sev- 
eral highly excited X-ray zones combined with less ex- 
cited photon-dominated regions are thought to be respon- 
sible for this behaviour. This combination yields Tk ~ 70- 
170 K and n(H2)= 10^"^ cm"^ globally, without indicat- 
ing a specific CO column density. Similarly to the case of 
M82, only the lower portions of these ranges are relevant 
here, i.e. those susceptible to reproduce the low-J CO tran- 
sitions. Another galax;y representing well AGN-dominated 
star formation activity, a lthough perhaps le ss exotic than 
Markarian231, is Arp220. ICreve et al.l (|2009l '). without pro- 
viding the CO SLED, show that the kinematic tempera- 
ture ranges from 40 to 60 K and the molecular gas vol- 
ume density n(H2)« 3 x 10^ cm~^, again without indicating 
a CO column density. Finall y, for the center o f the Milky 
Way, a normal spiral galaxy, iBavet et al. I (I2OO4') and refer- 
ences the rein showed th a t the CO SLED peaks at ^^C0(4- 
3), while ISanders et all (|l992l ) report Tk = 15-20 K and 
w(H2)= 10^~* cm~^, without indicating a CO column den- 
sity. 

Here, despite the lack of high-J observations to confirm 
the shape of our best-fit SLEDs, but assuming that the pre- 
dicted SLEDs shown in Figure 3] (black solid lines) are reli- 
able (see Section [4. 3p . we see that the predicted CO SLEDs 
of most of our ETGs (13/18) peak around the (3-2) or (4-3) 
transition. Taking into account the best-fit model parame- 
ters obtained, it is thus likely that most of our molecular 
gas-rich ETGs have similar gas excitation conditions, and 
have a star-formation activity similar to that present in the 
center of the Milky Way. IC0676 and NGC3607 show CO 
SLED turnovers located around the CO (5-4) transition 
that, when also taking into account the best-fit model pa- 
rameters (Tk, n(H2) and A'^(CO)) for those sources, suggest 
a closer proximity to Maffei2's characteristics. NGC1222 
and NGC6014 have molecular gas and CO SLED properties 
closer to those of NGC253, while NGC7465 seems in agree- 
ment with Markarian231. Again, however, although these 
reults are suggestive and interesting, we must be careful to 
not over-interpret them as they are extrapolations from rel- 
atively low-J CO lines only. With additional higher-J CO 
transition data, more precise comparisons and confirmation 
of these results will be possible. 



6 CONCLUSIONS 

Exploiting recently published molecular line ratios, we have 
modelled for the first time the CO and HCN-HCO"^ gas com- 
ponents in a sizeable sample of gas-rich early-type galax;- 
ies (18), this using a non-LTE theoretical method (i.e. an 
LVG code). We have separated these two gas components 
as the molecules studied are known to probe different gas 
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phases, hence a priori should trace different physical con- 
ditions. While the HCN and HCO^ lines alone leave the 
modelling of their corresponding gas phase degenerate, we 
are able to properly constrain the gas kinetic temperature, 
If2 volume density and CO column density of the CO gas 
component using the ^^CO(1-0, 2-1, 3-2) and ^^CO(1-0, 2-1) 
transitions. 

Most galaxies in our sample have a relatively cool 
(10-20 K) gas, a common H2 volume density (n(H2)~ 
10^~* cm~'^) and an average CO column density (Ai'(CO)> 
10^* cm~^). For the first time, we also provide the pre- 
dicted CO spectral line energy distributions of these 18 
ETGs, including the predicted SLED peak position. Taking 
into account the gas excitation conditions, derived physi- 
cal properties, and the model assumptions and limitations, 
most of the molecular gas in our gas-rich ETGs appears 
to host star-formation activity at a level similar to that 
present in the center of the Milky Way. Since their CO 
SLEDs show higher excitation conditions, a few galaxies 
(IC0676, NGC1222, NGC3607, NGC6014 and NGC7465) 
may have slightly more active star formation or additional 
sources of excitation (AGN, cosmic ray-dominated regions, 
etc). Strong 24, 60 and 160 ^m Spitzer fiuxes in NGC3607 
(compared to average values for ETGs) seem t o confirm this 
hypothesis l|Temi. Brighenti fc Mathews! [2OO9I '). 

The data used to constrain our models are averaged 
over several kiloparsecs and often over the entire galaxy. 
Local (sub-kpc) gas properties could thus be significantly 
difi^erent from the estimates presented here, and the values 
listed in Table [3] and [4] should be used with caution. Never- 
theless, these values are the best estimates of the molecular 
gas properties of ETGs available, and will remain so un- 
til spatially-resolved interferometric data of multiple species 
and transitions become available. 

This is also the first time that the physical properties 
are derived for two separate cold gas components in ETGs. 
This is a huge step forward in our understanding of current 
(residual) star formation in these objects. Despite the large 
uncertainties in the results derived for the HCN and HCO'*' 
gas component, our study nevertheless clearly shows that 
the physical properties of this component differ from those 
of the CO gas component, even in ETGs. This confirms the 
need for a separate treatment of these two gas components, 
with possibly different evolutionary paths. 

It is also essential to identify additional transitions of 
CO, HCN and HCO"*" likely to improve these results and 
better constrain the physical properties of the dense gas 
in ETGs in the future. We have identified here the mid-J 
CO transitions (e.g. ^^CO(6-5)) as particularly useful, with 
much constraining power. 
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APPENDIX A: RESULTS FOR THE HCN AND 
HCO+ GAS COMPONENT 

We present here the modeling results for the HCN and 
HCO^ gas component, for the galaxies where both lines are 
detected. Figure lAll shows the Ax'^ contours of the HCN- 
HCO^ models, similarly to Figure [2] for the CO gas com- 
ponent. Table lAll lists the best-fit model parameters. The 
HCN-HCO'*' results show greater uncertainties as the mod- 
els are underconstrained (one line ratio but three model pa- 
rameters). They should thus be considered indicative only 
and should be used with caution. 

This paper has been typeset from a T^j]X/ ffl^jX file prepared 
by the author. 
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Figure Al. Best-fit models for tlie HCN and HCO"*" gas component. Each panel shows the Ax^ 
HCO+ gas component as a function of the kinetic temperature Tk and H2 volume density n(H2), for three values of the CO column 
density 7V(CO) (hence three plots) centered on the best-fit value and indicated at the bottom of each plot. The model grid is shown with 
dark-grey dots (red dots in the online version), whereas the best-fit model of each galaxy (x^in "-"^ = 0) is labeled with a light-grey 
filled circle (green filled circle in the online version). The Ax^ contours and greyscales are typically for 1 to Scr confidence levels (i.e. 
Ax^ = 1-0, 2.7, 3.8, 5.0 and 6.6). The low a, high confidence level contours containing 'good' models are represented by the darkest 
areas. White zones correspond to models where there is no solution for fitting the single li ne ratio HCN(1-0)/HCO +(1-0). These models 
produce either negative opacities or negative integrated line intensities (see limitations in Ivan der Tak et"aDl2007l ). The galaxy name is 
indicated in the top-right corner of each panel. 



